Identical temperature dependence of the time scales of several 
linear-response functions of two glass-forming liquids 
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The frequency-dependent dielectric constant, shear and adiabatic bulk moduli, longitudinal thermal expan- 
sion coefficient, and longitudinal specific heat have been measured for two van der Waals glass-forming liquids, 
tetramethyl-tetraphenyl-trisiloxane (DC704) and 5-polyphenyl-4-ether. Within the experimental uncertain- 
ties the loss-peak frequencies of the measured response functions have identical temperature dependence over 
a range of temperatures, for which the Maxwell relaxation time varies more than nine orders of magnitude. 
The time scales are ordered from fastest to slowest as follows: Shear modulus, adiabatic bulk modulus, 
dielectric constant, longitudinal thermal expansion coefficient, longitudinal specific heat. The ordering is dis- 
cussed in light of the recent conjecture that van der Waals liquids are strongly correlating, i.e., approximate 
single-parameter liquids. 

PACS numbers: 64.70.P- 

Keywords: supercooled liquids; timescales; response functions; decoupling; ultra-viscous liquids 



A liquid has several characteristic timea^®. One is the 
Maxwell relaxation time determining how fast stress re- 
laxes tm = tj/Goc, where r\ is the shear viscosity and Goo 
the instantaneous shear modulus^. Other characteristic 
times arc identified by writing D = a 2 /td, in which D 
may be particle, heat, or the transverse momentum dif- 
fusion constant, and a is of order the intermolecular dis- 
tance. Further characteristic times are the inverse loss- 
peak frequencies (i.e., frequencies of maximum imaginary 
part) of different complex frequency-dependent linear- 
response functional For low-viscosity liquids like am- 
bient water the characteristic times are all of the same 
order of magnitude, in the picosecond range, and only 
weakly dependent on temperature. 

Supercooling a liquid increases dramatically its 
viscositjJsHUl most characteristic times likewise increase 
dramatically. The metastable equilibrium liquid can be 
cooled until the relaxation times become 10 10 -10 15 times 
larger than for the low-viscosity liquid, at which point 
the system falls out of metastable equilibrium and forms 
a glass. At typical laboratory cooling rates (K/min) 
the glass transition takes place when tm is of order 100 
seconda^S] 

Even though most characteristic times increase dra- 
matically when the liquid is cooled, they are generally not 
identical. Different measured quantities and different def- 
initions of the characteristic time scale lead to somewhat 
different characteristic times. A trivial example is the 
difference between the time scales of the bulk modulus 
and its inverse in the frequency domain, the bulk com- 
pressibility. More interestingly, some time scales might 
have quite different temperature dependence; this is often 
referred to as a decoupling of the corresponding micro- 
scopic processes. 
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Several relaxation time decouplings have been reported 
in the ultraviscous liquid state above the glass transition. 
Significant decoupling takes place for some glass-forming 
molten salts like "CKN" (a 60/40 mixture of Ca(N0 3 ) 2 
and KNO3 — ) , where the conductivity relaxation time at 
the glas s tra nsition is roughly 10 4 — 10 5 times smaller 
than T]v J 14 * 15 *- This reflects a decoupling of the molecu- 
lar motions, with the cations diffusing much faster than 
the nitrate ions^. A more recent discovery is the de- 
coupling of translational and rotational motion in most 
molecular liquids, for which one often finds that molecu- 
lar rotations are 10-100 times slower than expected from 
the diffusion time^KS. This is generally believed to re- 
flect dynamic heterogeneity of glass-forming liquidJ^HUl 
Angell in 1991 suggested a scenario consisting of "a series 
of decouplings which occurs on decreasing temperature" , 
sort of a hierarchy. He cautiously added, though, that 
"more data are urgently needed to decide if this repre- 
sents the general case'^. 

Some linear-response functions like the thermovis- 
coelastic and shear-mechanical ones are difficult to mea- 
sure reliably for ultraviscous liquidiP^H. To the best of 
our knowledge there are no studies of their possible de- 
coupling. This paper presents such data, together with 
conventional dielectric data. The purpose is to establish 
the order of relaxation times among the different quanti- 
ties and, in particular, to investigate whether or not they 
show a decoupling upon approaching the glass-transition 
temperature. 



I. EXPERIMENTAL RESULTS 

We have measured the complex, frequency -depe ndent 
dielectric constant e(ui), s hear modulus G(u>] ' 20 ' 23 \ adia- 
batic bulk modulus Ks(w'f^^, and longitudinal specific 
heat c/ (w}PH on two van der Waals bonded glass-forming 
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FIG. 1. Example of data from which loss-peak positions are 
found. Main part: Normalized imaginary parts of the follow- 
ing complex response functions: shear modulus G(uj), adi- 
abatic bulk modulus Ks(cj), dielectric constant e(uj), and 
longitudinal specific heat q(cj), as functions of frequency 
for 5-polyphenyl-4-ether (5PPE) at T — 252. 5K. Inset: 
Time-dependent longitudinal expansion coefficient ai(t) at 
T = 209. 2K for tetramethyl-tetraphenyl-trisiloxane (DC704). 



liquids tetramethyl-tetraphenyl-trisiloxane (DC704) and 
5-polyphenyl-4-ether (5PPE) (commercial vacuum-pump 
oils). For DC704 the time-dependent longitudinal ther- 
mal expansion coefficient ai (ipsa was a i so measured (see 
appendix [B] for details). We note that both liquids have 
linear-response functions that to a good ap proxima tion 
obey time-temperature superposition (TTS) 23 * 28 * 29 *, i.e., 
their loss-peak shapes are temperature independent in 
log-log plots. Moreover, the two liquids have only small 
beta relaxations, they rarely crystallize, and they are gen- 
erally very stable and reproducible — altogether these 
two liquids are very suitable for fundamental studies. 

Three of the measured quantities (a/, c;, and Kg) are 
closely related to one complete set of independent scalar 
thermoviscoelastic response functional: c p , a p , and kt 
(the relations are given in appendix [C]). Measurements 
of such a complete set of three scalar thermoviscoelastic 
response functions are rare, if at all existing for any glass- 
forming liquid. 

Figure [T] shows the normalized loss peaks as func- 
tions of frequency of G(oj), Ks(io), e(w), and ci{lo) for 
5PPE at 252. 5K. The inset shows one data set for the 
time-dependent a t (t) at T = 209. 2K for DC704. The 
frequency-domain data allow for direct determination of 
the loss-peak frequencies (i^ p ); the time-domain data 
were Laplace transformed to give an equivalent loss-peak 
frequency^. 

The measurements give both real and imaginary parts 
of the complex response functions, allowing us to calcu- 
late two other relevant characteristic frequencies, namely 
the inverse Maxwell time 1/(27ttm) (jj and G^ can be 
found from G(uj)) and the loss-peak frequency of the adi- 



abatic compressibility ks(uj) = l/Ks{w). Figure [2]( a) 
shows the temperature dependence of the seven charac- 
teristic frequencies for DC704 31 . The data cover more 
than nine decades in relaxation time (from T g and up). 
Clearly the time scales for the response functions follow 
each other closely. 

Figure [2](b) plots the characteristic frequencies in 
terms of a "time-scale index" defined as the logarith- 
mic distance to the dielectric loss-peak frequency^. For 
both liquids the time-scale indices are temperature in- 
dependent within the experimental uncertainty, that is, 
the characteristic time scales of the measured quantities 
change in the same way with temperature. 

This finding constitutes the main result of the present 
paper, showing that the time scales for theses response 
functions are strongly coupled, in contrast to the ob- 
served decoupling between transitional diffusion and 
rotation^HUl anc j a t variance with the scenario suggested 
by Angelffl 



II. DISCUSSION 

The fastest response function is the shear modulus. 
It has previously been reported that this quantity is 
faster than dielectric relaxation for a number of glass- 
forming liquids (see e.g. Ref. EH and Ref. [551 and refer- 
ences therein), a fact that the Gemant-DiMarzio-Bishop 
model explains qualitatively^. The dielectric relaxation 
is faster than the specific heat (this difference cannot be 
attributed to measuring q and not c j~) • For glycerol 
the same tendency has been report eoP^ZI, but with a 
fairly small difference in time scales of c p and e (« 0.2 
decades). However a consistent interpretation of dielec- 
tric hole-burning experiments on glycerol was arrived at 
by assuming that these two time scales coincide^. For 
propylene glycol the opposite trend has been reported^. 
Regarding volume and enthalpy relaxation there is like- 
wise no general trend in the literature; some glass-formers 
have slower enthalpy than volume relaxation, others the 
opposite^HSl. Clearly more work is needed to identify 
any possible general trends. However, such comparisons 
are difficult, as it requires a precision in absolute temper- 
ature at least better than IK between the experiments. 
This is very difficult to obtain, and it could be speculated 
that some of the contradicting results could be explained 
this way. An advantage of our methods is that the same 
cryostat can be used for all the measurements ensuring 
same absolute temperature (see appendix |B|) . 

What does theory have to say about the decoupling 
among relaxation functions and why some are faster than 
others? As mentioned, there are three independent scalar 
thermoviscoelastic response functions. There is no a pri- 
ori reason these should have even comparable loss-peak 
frequencies. Moreover, both the dielectric constant and 
the shear modulus are linear-response functions that do 
not belong to the class of scalar thermoviscoelastic re- 
sponse functions; these two functions could in principle 
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FIG. 2. (a) Seven characteristic frequencies (based on five 
measured quantities) as functions of temperature for DC704. 
The full curve is a fit of the dielectric data to a Parabolic 
function^ (see appendix [d] for details on this function), used 
for extrapolation of the dielectric data to low temperatures. 
Filled symbols (tm and ks) indicate that the quantity in ques- 
tion is not measured directly, but derived from one of the 
measured quantities. Error bars on a data were estimated by 
assuming an additive influence from the underlying spurious 
frequency dependence of the raw datcP, varying its influence 
on the loss peak from negligible to maximal. Equivalent data 
are given for 5PPE in appendix [A] 

(b) Time-scale index of all measured response functions (sym- 
bols as in the top figure) with respect to the dielectric con- 
stant for the two liquids DC704 and 5PPE. The time-scale 
index is plotted as functions of the dielectric loss-peak fre- 
quency (common X-axis), which represents the temperature 
(also given for each liquid). For both liquids the time-scale 
index for all quantities are temperature independent within 
the experimental uncertainty, that is, the measured quantities 
have the same temperature dependence of their characteristic 
time scales. 



have relaxation times entirely unrelated to those of the 
scalar thermoviscoelastic response functions. All in all, 
general theory does not explain our findings. 

As mentioned earlier, the two investigated liquids obey 
TTS to a good approximation and have very small (if 
existing) beta relaxations. In an earlier work (Ref. [29]) 
some of us noticed that the time-scale index between 
shear mechanical and dielectric relaxation is only signif- 
icantly temperature dependent for systems with a sig- 
nificant beta relaxation. In Ref. [42 it was shown for a 
number of systems (including systems with a pronounced 
beta relaxation) that the dielectric relaxation and the ag- 
ing rate after a temperature jump follow the same "in- 
ner clock"; these results were obtained at temperatures 
where the alpha and beta relaxations are well separated. 
Based on this, one might speculate that some (or all) of 
the temperature dependencies of the time-scale index ob- 
served in the literature could be due to the difference in 
the influence from the beta relaxation between the mea- 
sured response functions (see e.g. Ref. 03] for a compar- 
ison of the influence in shear mechanical and dielectric 
relaxation) . 



A. Comparison to a "perfectly correlating liquid" 

The class of so-called "strongly correlating" liquids was 
recently identified^. This class includes most or all van 
der Waals and metallic liquids, but not the covalently- 
bonded, hydrogen-bonded, or ionic liquids. In com- 
puter simulations strongly correlating liquids are charac- 
terized by strong correlations between constant-volume 
equilibrium fluctuations of virial and potential energy^. 
These liquids are approximate single-parameter liquids, 
i.e., they do not have three independent isotropic scalar 
thermoviscoelastic response functions, but to a good ap- 
proximation merely one^. A perfect single-parameter 
liquid obeys^ 



T a"(uj) k. 

iTp = 



e£(w) 



(1) 



where 7t p is a constant. Experimental evidence that 
DC704 is a strongly correlating liquid, was very recently 
presented in Ref. |47l 

We tested how well our results conform to the predic- 
tions for a perfectly correlating liquid. If the thermovis- 
coelastic response functions obey Eq. ([!]), the loss-peak 
frequencies should be identical at all temperatures. A 
single-parameter model liquid was constructed by assum- 
ing Eq. ([!]) to hold, with high-frequency limits and relax- 
ation strengths of the quantities chosen to be as close 
as possible to those measured for DC704. The quanti- 
ties ai(uj), cz(w), and ks(w) were calculated in the model 
by introducing shear-modulus data (see appendix [C] for 
details). 

Figure [3] compares the characteristic frequencies of q , 
ct\ , and Kg of our experiments to those of the model liq- 
uid. The order of the time scales for the model liquid 
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FIG. 3. Time-scale indices of cj, aj, and Kg for DC704 (blue 
symbols, left y-axis) and for the single-parameter model-liquid 
(red symbols, right y-axis) (see appendix [c] for details on 
the model calculation). The two axis have been shifted with 
respect to each other to give an overall best overlap of the 
indices. For the experimental data the time-scale index is 
given relative to the dielectric loss-peak frequency (i^ p , e ). The 
reported quantity is the mean over the index in Fig.[2jb), the 
error is based on the standard deviation (in the case of c; a 
mean over maximum and minimum based on the error bars on 
the loss-peak frequency was used). For the single-parameter 
model-liquid the index is with respect to the common loss- 
peak frequency of a p , c p , and kj- (fip,/)- The model liquid 
has the same order of the time scales as the real liquid, but 
the differences in the time scales are underestimated in the 
model. 



matches the observations, i.e., c; slower that cti and k$ 
faster than a;. However, the magnitude of the time-scale 
differences is significantly underestimated. 



III. SUMMARY 

We measured several complex frequency-dependent 
linear-response functions on the two van der Waals liq- 
uids DC704 and 5PPE. Within the experimental uncer- 
tainties the time scales of the response functions have the 
same temperature dependence, that is, the time-scale in- 
dices are temperature independent. The time scales are 
for both liquids ordered from fastest to slowest as fol- 
lows: Shear modulus, adiabatic bulk modulus, dielectric 
constant, longitudinal thermal expansion coefficient, lon- 
gitudinal specific heat. 

General theory does not explain why the time scales 
from some response functions couple very closely to each 
other, as is the case for the investigated response func- 
tions, and why others show a decoupling, as is observed, 
e.g. between transitional diffusion and rotatiorf^HIID The 
ordering of the longitudinal thermal expansion coeffi- 
cient, the longitudinal specific heat and the adiabatic 
compressibility can be rationalized by assuming that the 
liquids are stron gly co rrelating, i.e., approximate single- 
parameter liquids^ES in which certain sets of isotropic 



scalar thermoviscoelastic response functions have identi- 
cal time scales. 

More work is indeed needed to establish if the findings 
are general for van der Waals liquids, and to understand 
which response functions have temperature independent 
time scale index and which decouple. 
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FIG. 4. The six characteristic frequencies (based on the four 
measured quantities) as functions of temperature for 5PPE. 
Filled symbols (tm and Kg) indicate a quantity that is not 
directly measured, but derived from one of the other measured 
quantities. 



Appendix A: 5PPE characteristic frequencies 

Figure [4] shows the characteristic frequencies as func- 
tions of temperature^ for 5-polyphenyl-4-ether (5PPE) 
(the figure is equivalent to Fig. 2(a) in the main paper). 
Note that ai data are not available for 5PPE due to its 
larger dielectric relaxation strength^. 



Appendix B: Materials and methods 

The two substances investigated are both commercial 
vacuum-pump oils and were used as received without 
further purification. DC704 (tetramethyl-tetraphenyl- 
trisiloxane) was required from Aldrich as Dow Corn- 
ing ® silicone diffusion pump fluid [445975]. 5PPE 
(5-polyphenyl-4-ether) was required as Santovac ® 5 
Poly phenyl ether. 

With the exception of a\ (t) all quantities were mea- 
sured in the same cryostatF^GS ensuring the same abso- 
lute temperature. The thermal expansion coefficient data 
were taken in a different cryostat, for which the absolute 
temperature was calibrated using the liquid's dielectric 
relaxation time in a temperature range where data exists 
from both cryostats. ft was not possible to measure ai 
for 5PPE, because on was measured by a method that 
requires a very small dielectric relaxation strength^. 



Appendix C: The perfectly correlating (single-parameter) 
model liquid 

As stated in the main text a strongly correlating 
liquid^ is characterized in computer simulations by 
strong correlations between the constant-volume equilib- 
rium fluctuations of virial and potential energy (correla- 
tion coefficient above 0.9). Such a liquid is an approxi- 
mate single-parameter liquid, i.e., it does not have three 
independent isotropic thermoviscoelastic response func- 
tions, but merely one (to a good approximation) 46 . 

The purpose of the model calculation is to investigate 
to which extent our findings regarding the difference in 
characteristic frequency between the thermoviscoelastic 
response functions c/, ai, and Kg are consistent with 
the hypothesis that the investigated liquid tetramethyl- 
tctraphcnyl-trisiloxanc (DC704) is a strongly correlating 
liquid. For this purpose a perfect single-parameter model 
liquid is constructed and the loss-peak frequencies of the 
measured response functions are calculated. 

At all frequencies a perfect single-parameter liquid 
obeysP 



: 7Tp 



;( W ) 



(CI) 



where "fTp is a constant. By the Kramers-Kronig rela- 
tions Eq. (CI I implies that the relaxation strengths are 
related by 

T Aa n 



Ac, 



- = iTp = -T • 

Aa v 



(C2) 



It follows that the full response functions can be writ- 
ten as 



a p (uj) = Aa p f(uj) + a poo 
T 

c p {uj) = Aa p f(uj) + c pa 

iTp 

k t (uj) = j Tp Aa p f(uj) + K To 



(C3a) 
(C3b) 
(C3c) 



where /(w) is a normalized complex susceptibility, defin- 
ing the common time scale as well as the shape of the 
relaxation functions. a poo , c poo , and ktoo are the high- 
frequency limits of the response functions, and Aa p is 
the relaxation strength of a p . 

The measu red qu antiti es q, ai, and us ar e related to 
those of Eq. ((Cl b^^. 



ai(u>) 

Cy(w) 

ci(u)) 



a p (ui) 



1 

c p (u>) - 

Cp(uj) 



4/3 G(lj)k t (uj) 
_ T a p (uj) 2 

Kt(oj) 



+4/3 G{uj) 
\/kt{u) +4/3 G{uj) 



Cy(w). 



(C4a) 
(C4b) 
(C4c) 
(C4d) 



Thus if the shear modulus is known, it is possible to 
calculate the measured quantities based on the quantities 
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which are related in the single-parameter liquid: a p , c p , 
and kt- 

A single-parameter model liquid with properties close 
to DC704 was constructed in the following way (the pa- 
rameters used are listed in Table |l| : 

• The high-frequency limits and relaxation strengths 
of dp, c p , and kt were calculated based on the high- 
and low- frequency limits of of q, Ks, on, and G as 
given in Ref. H7] This can be done, because Eq. 
( C4 ) can be inverted analytically. 



• 7Tp was estimated from Eq. (C2). DC704 is not a 
perfectly correlating liquid, so the right- and left- 
hand sides do not give exactly the same value. They 
differ by ps 10%, which is below their respectively 
uncertainties, and the average was used. 

• The normalized complex susceptibility f(u>) was 
chosen to follow a "generalized BEL model" , which 
we found to des cribe th e alpha relaxation in many 
liquids very wel P 4 ' 43 ' 5 " !; 



i 



(C5) 



-k(iur) 



The parameters were estimated by fitting a modu- 
lus version of the model to the shear modulus data. 

• The shear modulus was assumed to follow a mod- 
ulus version of Eq. (C5), having G^ = 3/5AKs 
(in Ref. [23] this relation was shown to apply for 
DC704), and a time scale that is one decade faster 
than f(ui) (as seen in Fig. 2 of the main paper). 

• 10 5 realizations of the model liquid were calculated 
with the parameters 7t p , Aa p , a poo , c poo , and ktoo 
chosen independently from Gaussian distributions 
with mean given as in Table [I] and a standard devi- 
ation of 20%. This was done in order to investigate 



TABLE I. Parameters for the perfect single-parameter liquid 
mimicking the properties of DC704 at the temperature 214K. 
The values are based on the data presented in Refs. [23] and 

H7J 



7Tp 

A a p 

C^poo 



Goo 

logio 



2.5 
3.2 
1.4 
1.4 
1.9 



10 _7 m 3 KJ 
10" 

io- 

10 6 J(m 3 K) 



- 4 1/K 
- 4 1/K 



3/5AK s 
1.1 



a 
k 



1 + 
0.44 
1.05 



(iuT)- 1 +k(iur)- 



the robustness of the results with respect to the in- 
fluence of the absolute levels. The chosen standard 
deviation is larger than the estimated uncertainties 
on the quantities, ensuring that the results repre- 
sent a "worst case scenario". For each realization 
the loss-peak positions of q, ai, and k$ (and addi- 
tionally cy) were calculated. 

Figure [5] shows the time-scale index between the loss- 
peak frequencies of the calculated quantities with respect 
to the loss-peak frequency of f(u>) (Eq. (C5|). That is, 



Fig. [5] shows the loss-peak positions of c/ , cy , on and n$ 
when the underlying liquid is a perfect single-parameter 
liquid in which c p , a p , and kt have the same loss-peak 
position. It can be seen that the influence on the loss- 
peak position of measuring the longitudinal specific heat 
and thermal expansion coefficient, instead of their iso- 
baric counterparts, is very small. Compressibility, on the 
other hand, is affected more by measuring the adiabatic 
version instead of the isothermal. 



Appendix D: Choice of extrapolation function 

In order to extrapolate the dielectric data to the tem- 
peratures of the thermal expansion coefficient, a fitting 
function is needed. Inspired by the work of Elmatad et 
alW^ five fitting functions were investigated (Table [ll| : 
The Vogel-Fulche r-Ta mmann (VFT) functiorP^^, the 
Avramov functiorP^IHIE the Parabolic function^, and 
two double-exponential functions: The MYEGA function 
suggested by Mauro et atW^ and the similar FF2 function 
of Hecksher et alW^. All five functions have three fitting 
parameters, one of which is the microscopic attack fre- 
quency (i/ ). 

Each function was fitted to the loss-peak position from 
the dielectric constant, using a least-square minimiza- 
tion. Figure [6] illustrates the quality of the fits, which for 
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FIG. 5. Time-scale index for cj, oti, ks, and cv relative to the 
loss peak position of the common function /(u>) entering Eq. 

del. 



TABLE II. Fitting functions. The table presents the five functions (see main text for references) which were tested for the 
ability to extrapolate the dielectric loss peak to the temperatures where the thermal expansion coefficient was measured. The 
parameters reported derive from fitting the full data set. The "Parameter change" is the relative change in the parameters 
between fitting to the 12 highest temperatures (down to 224K) and to the full set of 24 points (down to 211. 5K). 

v\ p {T) Parameters Parameter chage 

VFT i/o exp (-jSjt) A = 3510K, T = 148.9K, log 10 (i> ) = 22.1 A : 62%, T a : -10%, log 10 (z/ ) : 22% 



T-T 0/ 

§)") n = 5.0, B = 437.0K,log 10 (!/ ) = 14.0 n : —19%, B : 17%, log 10 (i/ ) : 17% 



Parabolic v exp (- J 2 (± - J = 3740K 2 , T = 299.6K, log 10 (i/ ) = 9.5 J : -8%, T : 5%, log w (u ) : 11% 

Double exp (MYEGA) v exp (-f exp (§)) K = 294.4K,C* = 722.5K, log 10 (i/ ) = 16.2 K : 344%, C : -28%, Iog 10 (i/ ) : 22% 
Double exp (FF2) ko exp (-Aexp (^)) ,4 = 0.650, T 2 = 893.5K, log 10 (i/ ) = 17.0 A : 393%, T 2 : -25%, log 10 (v ) : 25% 
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log 10 (vip, F it/Hz) 

FIG. 6. Quality of five fitting functions of the loss-peak fre- 
quency as functions of temperature. The main figure shows 
the result of evaluating fits of the five different function to 
the full data set (see Table [II] for details on expressions and 
parameters and the main text for references) , plotted against 
the measured loss-peak frequencies (for clarity the plots have 
been offset by an integer constant on the y-axis). For a per- 
fect fit the points fall on the shown black lines. The (barely 
visible) dashed green lines show the ability to extrapolate; it 
is the result of evaluating a fit including only data down to 
log 10 (^i p ) = 0.67 (corresponding to 220K). The insert shows 
the best fit attack frequencies, i/o, as functions of the lowest 
temperature included in the fit (see the main text for details). 
This figures shows that all five functions fit data well, but the 
Parabolic function is the best function to used for extrapola- 
tion of data to lower temperatures. 



all five functions are excellent. The fitting parameters 
are given in Table [IT] In order to explore to which extent 
the functions are useful for extrapolating, the stability 
of the parameters were investigated when adding points 
to the fit. All functions were fitted to the 12 points at 
the highest temperatures. Points were then added to the 
data set one by one and the functions re-fitted. The rel- 
ative change in the fitting parameters from the initial to 
the final data set is reported in Table \H\ and the insert 
in Fig. [6] shows the fits for the attack frequencies (vq). 
A clear convergence of the parameters was not observed 
for any of the functions, but the Parabolic function has 
the most stable parameters (however, the Avramov func- 
tion is only slightly worse and it has a more physically 
realistic attack frequency) . Finally, the ability to extrap- 
olate was investigated by fitting to the first 16 points and 
then evaluating the fitted function at the lower temper- 
atures, as shown by the dashed green lines in Fig. [6j). 
The extrapolation spans 2.5 decade in loss-peak fre- 
quency, comparable to the extrapolation needed to get 
to the lowest temperature of the thermal expansion co- 
efficient. Again the Parabolic function is best, but not 
perfect. Based on this the Parabolic function was chosen 
as extrapolation function. 



